Spin-i Heisenberg-Antiferromagnet on the Kagome Lattice: High 
Temperature Expansion and Exact Diagonalisation Studies 

N. Elstner and A. P. Young 
Physics Department, University of California, Santa Cruz, CA 95064 

Abstract 

For the spin-^ Heisenberg antiferromagnet on the Kagome lattice we cal- 
culate the high temperature series for the specific heat and the structure 
factor. A comparison of the series with exact diagonalisation studies shows 
that the specific heat has further structure at lower temperature in addition to 
a high temperature peak at T ~ 2/3. At T = 0.25 the structure factor agrees 
quite well with results for the ground state of a finite cluster with 36 sites. 
At this temperature the structure factor is less than two times its T = oo 
value and depends only weakly on the wavevector q, indicating the absence of 
magnetic order and a correlation length of less than one lattice spacing. The 
uniform susceptibility has a maximum at T ~ 1/6 and vanishes exponentially 

for lower temperatures. 
PACS numbers: 75.10.Jm 
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I. INTRODUCTION 



For a long a time it has been speculated that low-dimensional quantum spin antiferro- 
magnets may have magnetically disordered ground states. Most attention has been focused 
on the spin- 1 Heisenberg antiferromagnet (HAFM) on the square lattice due to the close re- 
lation of this model with the problem of high temperature superconductivity. However, it is 
now well established that the HAFM on this particular lattice has an ordered groundstateS. 
The first system for which a disordered groundstate was proposed is the HAFM on the 
triangular lattice!. In recent years this model has been the subject of intensive numerical 
investigations. Although most results indicate that the system remains ordered at T = 0, it 
seems that the sublattice magnetisation and the spin stiffness are significantly smaller than 
for the square latticei~i 

So far the best candidate for a magnetically disordered system is the Kagome structure: 
a triangular lattice with a triangular basis. The vectors of the underlying triangular Bravais 
lattice are 

ei = 2(1,0) , e 2 = (l,v / 3) (1) 

while basis vectors indicating the coordinates of the three sites in the triangular unit cell 
are 

b x = (0,0) , b 2 = (l,0) , b 3 =^,^ . (2) 

The units are chosen so that the nearest neighbour distance equals unity. The structure is 
shown in fig.l. Because the Kagome lattice is not a Bravais lattice, but has three sites per 
unit cell, the structure factor is a 3 x 3 matrix. 

Classical spins on the Kagome structure are frustrated, as they are on the triangular 
lattice, but the coordination number is smaller (4 for the Kagome structure instead of 6 for 
the triangular lattice). Even more important may be another difference: while the ground 
state on the triangular lattice is degenerate only with respect to global rotations in spin 
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space, the classical ground states on the Kagome lattice have a local degeneracy, which 
results in a finite ground state entropy. Fluctuations around magnetically ordered states 
have a dispersionless zero-energy modeB. 

Numerical studies using series expansions^ and exact diagonalisation techniques^ have 
convincingly shown the the ground state of the spin-| HAFM on the Kagome lattice has no 
long range magnetic order. The exact diagonalisation studies!'! find a very rapid decay of 
the spin-spin correlations indicating a correlation length £ of only about one lattice spacing 
and a finite spin gap A ph 0.25 J. There exist a number of proposals for a disordered ground 
state0~E!. Large-N expansions for the SU(iV) and Sp(iV/2) generalisations!^ of the HAFM 
predict a ground state with spin-Peierls order for SU(iV) or a spin liquid for Sp(iV/2). Up 
to now there exist no results from numerical calculations that confirm or contradict any of 
these proposals. 

Apart from being a theoretical toy model for a disordered quantum spin system there 
exists at least one possible realisation of the spin-| HAFM on the Kagome structure; the 
second layer of 3 He atoms absorbed onto graphite at a particular coverag It is, however, 
possible that this is too simple a model and a realistic description should include other spin- 
exchange interactions!!^. Experiments on this system found a peak in the specific lieatQ, but 
the total change in entropy per site between T = oo and an extrapolated value for T = 
accounts for only one half of ln(2) the expected valuefor a S=| system. This suggests that 
there is a large number of low-lying states, which could contribute to additional structure, 
such as a second peak, in the specific heat at very low temperatures. 

Exact diagonalisation of a 12-site cluster^ on the Kagome lattice and simulations us- 
ing the decoupled- cell Monte Carlo technique! found such a peak. Simulations of larger 
systems^! using the forced occillator method found only one high temperature peak and an 
almost linear T-dependance of the specific heat at temperatures below this single maximum. 
Based on this observation it was concluded that the double peak structure reported for the 
12-site system is due to finite size effects. 

Further insight in the behaviour of the system can be gained by calculating finite temper- 



ature properties. However, the most powerful technique, Quantum Monte Carlo simulations, 
breaks down due to the sign problem for frustrated spin systems. In this paper we, therefore, 
follow a different approach and present results from high temperature expansions as well as 
exact diagonalisation studies. 

Section 2 contains the results for the specfic heat, for which the series has been calculated 
up to 16th order in J jk^I '. We analyse the series using the method of Pade approximations 
and compare the results with data from exact diagonalizations studies. Based on entropy 
arguments we will show that the low temperature structure of the specific heat occuring in 
the finite cluster calculations^ cannot be a spurious finite size effect. In the third section 
we present results for the spin-spin correlation function, which has been calculated up to 
14th order in J/k^T. Extrapolations give quantitative results down to T »s 0.25 J/ k-Q. 
At this temperature the largest eigenvalue of the structure factor matrix S a) p(q) is nearly 
independent of the wavevector q and is less than a factor of two larger than at T = oo. 
Our values for the structure factor at this temperature are close to exact results for the 
groundstate of a finite cluster with 36 sites!. Our conclusions are summarized in section 4. 



II. SPECIFIC HEAT 

The Hamiltonian of the Heisenberg model is given by 

H = J S i • S i J > , (3) 

(ij) 

where the sum runs over all pairs of nearest neighbours on the Kagome lattice. We calculated 
the high temperature series using a linked cluster expansion^! up to order 16 in 1/T (from 
now on we will set the exchange coupling J and the Boltzmann factor k-Q equal to unity). 
The series coefficients are given in table 1. The series was extrapolated beyond its radius 
of convergence by the method of Pade approximants. For a power series F(x) we form the 
Pade approximants 

WM] = £m (4) 



where Pl{x) and Qm{x) are polynomials in x of order L and M respectively. The coefficients 
of the two polynomials are determined by the condition that the expansion of [L/M] has to 
agree with the series F(x) up to order 0(x L+M ). Asymptotically a [L/M]-Pade approximant 
has the behaviour 

Jim [L/M] oc x L ~ M . (5) 

In case of high temperature expansions x is the inverse temperature. Because the specific 
heat must vanish at T = we restrict the Pade analysis to approximants with M > L. 

A number of approximants were obtained this way and are plotted in fig. 2. The curves 
for different approximants remain consistent with each other down to T = 0. Furthermore, 
our findings agree with results from simulation studies on small clusters (of up to 18 sites)0. 
In particular there is only one peak at T « 2/3. We calculated the total change in entropy 

r°° C 

AS = l -dT (6) 

and found 

A 9 

=f«0.61n(2) (7) 

from integration of the Pade approximants for the specific heat. The [L/L] Pades for the 
entropy, which by construction go to a finite value at zero temperature, are in agreement 
with this result. Experiments on 3 He films0 absorbed on graphite reported a total change 
in entropy per site of only about ~ ln(2) which is very close to our results. However, this 
is not enough to show that this model is appropriate to describe these experiments. There 
exist other proposals^ that give similar findings. 

We do not, however, expect the large ground state degeneracy corresponding to a T = 
entropy per site of 0.4 In 2 implied by Eq. (7). Presumably tunneling removes the high 
ground state degeneracy which occurs in the classical model, leading to many low lying (but 
split) states for the quantum case. This would give additional structure to the specific heat 
at low temperatures where the high-T series do not converge and where the experiments 
have not yet been performed. 



We therefore also calculated the complete spectrum for finite systems with N = 12, 
15 and 18 sites using standard diagonalisation routinesEOL We chose two different clusters 
with 18 sites. The two possibilities, referred to as 18a and 18b respectively, are shown in 
fig. 1. Cluster 18b is the one used in the exact diagonalisation studies by Zeng and Elseii 
and the simulations performed by Fukamachi and NishimoriEl. Because translations are 
the only possible symmetry operations of this cluster we were not able able to reduce the 
dimension of the Hilbert space to a size that routines for complete diagonalisation could 
be used. However, we calculated a large number of the low-lying eigenvalues by a Lanczos 
algorithmic, which gives the accurate specific heat at the low temperatures of interest. Due 
to the larger symmetry of cluster 18a, we were able to calculate all eigenvalues of this 
particular system as well as for the N=12 and 15 clusters. 

Results for the specific heat are presented in fig. 3 together with one of the Pade ap- 
proximants for comparison. The two different methods, high temperature expansions in the 
thermodynamic limit and exact diagonalisation for finite clusters, are in excellent agreement 
with each other down to T « 0.3, i.e. below the high temperature peak. Thus, results from 
the finite cluster calculations are already in the thermodynamic limit for T > 0.3 

At lower temperatures, however, the two method give completely different results. The 
finite clusters with an even number of sites develop a sharp peak in the specific heat, while 
the 15 site system has at least a clear shoulder. In simulationsEl the shoulder for the 15 
site system appeared to be much less significant and no peak was observed in the 18 site 
system. This led to the conclusions that the low temperature peak reported earlier0 was a 
finite size effect. Here we see that there is additional structure in the specific heat at low 
temperatures for larger sizes. We expect that some structure persists in the thermodynamic 
limit leading to a vanishing entropy as T — > 0. The precise form of this structure, e.g. 
broad shoulder or second peak, however, is difficult to deduce from our results, because 
the finite size corrections show a large even-odd asymmetry. Exact diagonalisation studies 
of finite clusters! find that the lowest triplet (quadruplet) for finite systems with an even 
(odd) number of sites N has an excitation energy that remains finite in the thermodynamic 



limit resulting in a spin gap A pa 0.25. The low-lying states that give rise to the additional 
structure in the specific heat are singlets (doublets) for even (odd) N. 

III. STRUCTURE FACTOR AND UNIFORM SUSCEPTIBILITY 

In addition to the specific heat we also calculated the high temperature series for the 
spin-spin correlations. Only 14 terms were determined for this series, because many more 
clusters contribute to this expansion than for the specific heat. From the series for the 
correlations in real space we calculated the structure factor. As already mentioned in the 
introduction the Kagome structure is not a Bravais lattice so the structure factor is a 3 x 3 
matrix, given by 

S*A<1) = E expHq • (R + bg - b Q )] (S z {b a ) ■ S Z (R + b^) (8) 

R 

where R is summed over Bravais lattice vectors formed from e x and e 2 in eq. (1) and the 
b a , a — 1, 2 or 3 are given by eq. (2). Harris et al.S showed that up to sixth order in 1/T 
the largest eigenvalue of this matrix is independent of the wavevector q. This effect is due to 
the geometrical properties of the lattice. This degeneracy is broken in 7th order for quantum 
spins and in 8th order for the classical case. For certain values of the wavevector with high 
symmetry the eigenvectors are independent of T so the eigenvalues can be obtained as series 
in 1/T. The results for the largest eigenvalue at q = 0, q = 0) (corner of the Brillouin 
zone) and q = ^(0, 1) (center of an edge of the Brillouin zone) are given in table 1. It is 
remarkable that the values for the different wavevectors differ by only about 5% indicating 
only weak dispersion. The q-value with the largest coeffcient in the structure factor series 
changes at each order of the expansion, see table 1. This questions the conclusion in the 
classical limit@, for which a tendancy towards selection of Q = ^(1,0) was reported based 
on an eight term series. The next order of the expansion may already change that. 

For wavevectors not at high symmetry points in the Brillouin zone we first calculated 
the Pade approximants for all nine matrix elements of S at p(q), evaluated them at a fixed 
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temperature and then diagonalized the matrix. In fig. 4 we present results of a scan through 
the Brillouin zone at temperature T = 0.25. For this data we first made a transformation 
to the new variable u = tanh(//2T) where / = 1/8 and performed the Pade analysis in the 
series for this variable. We found that the transformed series behaved better then the original 
one. Figure 4 shows that the largest eigenvalue is nearly independant of q. Although the 
different [L/M] Pades have a weak dispersion, there is no clear tendancy towards selection 
of a particular q- value. We estimate that the largest eigenvalue of S a> p(q) is given by 

4S' ma:E (q,T = 0.25) = 1.72 ±0.04, (9) 

and any dispersion is smaller than the error bars. For comparison: the corresponding result 
at T = oo is AS max = 1. 

In order to check whether our results give a correct picture of the low temperature 
magnetic properties we calculated the structure factor for the groundstate of a finite cluster 
with 36 sites. The groundstate of this cluster has wavevector q = 0. The correlations in real 
space are given in ref. 12. For the 36 site system the allowed wavevectors are q = (0, 0) , 
|(1, 0) , ^r(l, 0) , ^(0, 1) and others related by symmetry to these. The results are plotted 
in figure 4. The agreement between the results obtained from the series analysis at T = 0.25 
and the ground state of the 36 site cluster is remarkable. It indicates that the results of the 
series analysis remain qualitatively unchanged down to T = 0, i.e. there is no divergence of 
the structure factor which would indicate magnetic order. The structure factor at T = is 
roughly a factor of two larger than at T = oo and has only moderate dispersion, indicating 
a correlation length of less than one lattice spacing. 

Results for the uniform susceptibility are plotted in figure 5 using the series for the sus- 
ceptibility in table 1. We show only results for finite cluster with an even number of sites 
N. As already mentioned at the end of the previous section the lowest states for a finite 
system with an odd number of sites are doublets which results in a strong finite effect for 
the susceptibility: lim^o X — ^ esu ^ s from the high temperature expansion and from 

finite cluster calculations again agree down to T ^ 0.3. Above this temperature the suscep- 
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tibility is significantly smaller than the Curie- Weiss susceptibility. At low temperatures the 
susceptibility appears to vanish exponentially, because the low-lying states are all singlets. 
As mentioned earlier the spin gap A is estimated^ to be A « 0.25. Because the the S z = ±1 
components of the lowest triplet give identical contributions to the susceptibility one expects 
the maximum in x to occur at temperature T max »i A/(l + ln(2)). This can be verified for 
the finite clusters with N sites by taking the finite size value A(JV). The maximum for the 
susceptibility is given by Xmax = 0.14 — 0.15 with only a rather weak iV dependance. 

Experiments on 3 He films absorbed on graphite find a cusp in the susceptibility near 
1 mKEl. Unfortunately our data do not allow to give quantitative estimates for the position 
and peak value of the susceptibility at such low temperatures, because the series expansion 
is no longer reliable and the finite cluster calculation suffer from finite size effects. It is very 
likely that in this temperature regime, T <C 1, the spin-| HAFM is no longer an appropriate 
model for 3 He films and additional interactions have to be taken into accountEi. 



IV. CONCLUSIONS 

We have presented results from high temperature expansions and exact diagonalisation 
studies for the specific heat and the structure factor of the Heisenberg antiferromagnet on 
the Kagome lattice. Our main result is that the specific heat has additional structure a 
second peak or possibly a shoulder at very low temperatures in addition to a peak at higher 
temperature (T rs 2/3). This had been conjectured earlier0 but subsequent simulationslll 
had seemed to contradict it. Our conclusions that this unusual low temperature behaviour 
exists comes from the follwing observations: high temperature series expansions, which do 
not find a second peak, obviously cannot account for a significant amount of entropy. At 
temperatures T < 0.25 finite cluster calculations show much more structures (a peak for even 
number of sites N, a significant shoulder for N odd) in the specific heat. For T > 1/4 the 
results of both methods agree, showing that in this regime the finite cluster calculations give 
the correct results in the thermodynamic limit. Therefore, the specific heat below T < 0.25 
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has to be much larger than previously reported to get the entropy right, i.e. zero as T — > 0. 
Thus, what is seen in the finite cluster calculations is not a spurious finite size effect. It 
would be interesting to measure the specific heat at somewhat lower temperature to see if 
additional structure appears. 

The series expansions for the structure factor at T = 0.25 give results that are very 
similar to the results for the ground state in exact diagonalisation studies!. Both indicate 
that correlations fall of rapidly with distance. The structure factor is less than a factor of 
two larger the infinite temperature value even at T = 0.25 . Any possible dispersion is to 
weak too be clearly identified. All this indicates, that the Heisenberg antiferromagnet on 
the Kagome lattice has no long range magnetic order. 

For the susceptibility we find a maximum around T ~ I and an exponential drop at 
lower temperatures. 
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Table 1 



N 


AC 


4x 











1 





4 


2 


48 


-32 


3 





192 


4 


-9792 


-384 


5 





-1280 


6 


4106880 


-155136 


7 


-5193216 


2711184 


8 


-2927834112 


56705024 


9 


11470159872 


-1716811776 


10 


3193027983360 


-47711784960 


11 


-26121748561920 


2004747075584 


12 


-4944246830899200 


55843726884864 


13 


70892246893658112 


-3367208347123712 


14 


10284867640404983808 


-88720801213743104 


15 


-234226245436710912000 


7723917022263705600 


16 


-27538523697287477329920 





For each quantity A we define coefficients, a n , by A = J2n=o (f j • The table shows 
the values of the a n for the specific heat C and the uniform susceptibility x 
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Table 2 



N 


4S max {q = 0) 




4Wq=f (1,0)) 





1 


1 


1 


1 


2 


2 


2 


2 


4 


4 


4 


3 


-72 


-72 


-72 


4 


-448 


-448 


-448 


5 


11872 


11872 


11872 


6 


122368 


122368 


122368 


7 


-4503872 


-4494912 


-4493120 


8 


-61508608 


-61640704 


-61749760 


9 


3088187904 


3072426496 


3066299904 


10 


48686666752 


49120629760 


49365544960 


11 


-3348876193792 


-3321064779776 


-3305454742528 


12 


-54711166472192 


-55953620766720 


-56467096948736 


13 


5268689812606976 


5210400086315008 


5168484556144640 


14 


80271309635928064 


84156994384281600 


85379333782306816 



For each quantity A we define coefficients, a n , by A = J2n=o fff (f J • The table shows 
the values of the a n for the largest eigenvalue ^^(q) of the structure factor matrix at wave 
vectors q = (0, 0) , q = -^(0, 1) and q = f (1, 0) 
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FIGURES 

FIG. 1. The Kagome lattice with the finite clusters used in exact diagonalisation studies of the 
specific heat. 

FIG. 2. [L/M] Pade approximants for the specific heat C obtained from a 16 term high tem- 
perature expansion. 

FIG. 3. The specific heat C calculated from exact diagonalisation of finite clusters compared 
with results from the Pade analysis of the high temperature series. For cluster 18b only the 
low-lying eigenvalues have been determined and so the results are only presented for T < 0.4, the 
region where they are valid. 

FIG. 4. The eigenvalues of the structure factor matrix at T = 0.25. The Pade analysis was 
done in the new variable u = tanh(//2T) where / = 0.125. The asterix marks the eigenvalues of 
the structure factor for the ground state of a finite cluster with 36 sites. The data were obtained 
using the results for the spin-spin correlations of this cluster given in ref. 12 

FIG. 5. Uniform susceptibility % from finite clusters and a Pade analysis of the high temperature 
series. The Pade analysis was performed for Tln(Tx). The data shown were obtained from the 
[7/5] Pade approximant. The Curie- Weiss susceptibility 4xcw = — @) with = — 1 is shown 
for comparison. The first two terms in the high temperature series expansions for x an d Xcw are 
equal, so x — Xcw at sufficiently high T. 
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